Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 



(MN MfeK style file vl.4) 



Cluster abundance and large scale structure 

Jiun-Huei Proty Wu* 

Astronomy Department, University of California, Berkeley, 601 Campbell Hall, Berkeley, CA 94720-341 1, USA 



1 February 2008 



ABSTRACT 

We use the presently observed number density of large X-ray clusters and the linear 
mass power spectra from galaxy surveys to constrain the amplitude of matter density 
perturbations on the scale of 8/i _1 Mpc (as), and the redshift distortion parameter ((3), 
in both open cosmologies and flat models with a non-zero cosmological constant. The 
best fit to the observed mass power spectra gives n = 0.84 ± 0.67 and T = 0.271$'^, 
with the theoretically expected degeneracy V = 0.247Texp(1.4n) = 0.220to;o3i (all at 
95 per cent confidence level). These are consistent with the recent CMB results. Based 
on this, we then calculate the cluster-abundance- normalized as, using different models 
of mass function. The models considered are by Press & Schechter (PS; 1974), by Sheth 
& Tormen (ST; 1999), and by Lee & Shandarin (LS; 1999). The last two incorporate 
non-spherical gravitational collapse, and the as based on these two models are signifi- 
cantly lower. This lower normalization results from the larger mass function within the 
scale range of our interest. In particular, we combine the results of these two models 
to yield ct 8 (st+ls) = 0.477f2^ , where a = — 0.3 — 0.17£!^q 4 — 0.13f2Ao- In our analysis, 
we also derive the probability distribution function of cluster formation redshift using 
the Lacey-Cole formalism (1993), but with modifications to incorporate non-spherical 
collapse. The origins of uncertainties in our as results are also investigated separately, 
with the main contributer being the normalization in the virial mass-temperature re- 
lation. From the PSCz power spectrum alone and using T' = 0.220to'o3i as the prior, 
we also obtain for the IRAS galaxies a 8 (i) = 0.78 ± 0.06 (at 95 per cent confidence 
level) . By combining this with the as result, we are able to constrain the redshift dis- 
tortion parameter [3\, which is in turn lower in the non-spherical-collapse models. We 
found /3i(st+ls) — 0.613f2^Q 4 °- 16 ( n '>>°+ nA <>)_ This is more consistent with the recent 
observations than the result based on the PS formalism. 

Key words: galaxies: clusters - large-scale structure of Universe - X-rays - cosmol- 
ogy: theory 



1 INTRODUCTION 

One of the most important constraints on models of struc- 
ture formation is the observed abundance of galaxy clusters. 
Because they are the largest virialized objects in the uni- 
verse, their abundance can be simply predicted by linear per- 



turbation llicuiy. In lite lileiHlme, Lite cluster abuud 



spectrum of mass perturbations, the so-called 'mass func- 
tion' (the differential number density of collapsed objects as 
a function of mass) , the probability distribution function of 
cluster formation redshift, and the virial mass-temperature 
relatio n. Although this has been studied by many other au- 
thors flViana fc Liddle 1999|; |Eke, Cole fc Frenk 1996fc |Wang 



Steinhardt 1998| ; |Borgaiii ct a.l. I ')')']: |Picrpaoli. Scott & 



(White 
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White 2000), we revisit this problem with more caution to 



Viana fc LTddle 1999|; |Kitayama fc Suto 19971 ; [Wang fc Stein, 
hardt 199^ ; |Avclmo, Wu fc Shellard 2000J )' Here, we shall 
consider its constraint on the amplitude of matter density 
perturbations on the scale of 8/i _1 Mpc, as, and the red- 
shift distortion parameter, /3, in the standard inflationary 
models. This requires the knowledge about the linear power 
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incorporate the recent progress in measuring the linear mass 
power spectrum, as well as understanding the mass func- 
tion. Based on these new developments, we shall investigate 
as and f3 in both open cosmologies and flat models with a 
non-zero cosmological constant. 

The structure of the paper is as follows. In section ^[ 
we first investigate the observed linear mass power spectra. 
This includes the estimation of erg, the spectral index n, and 
the shape parameter F. They parameterize the linear mass 
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power spectrum predicted by the standard inflationary mod- 
els. Within such models, we reveal the degeneracy between n 
and F, and thus reparameterize the linear mass power spec- 
trum with only a single parameter, the degenerated shape 
parameter V' . The IRAS cr g m is also estimated here. In sec- 
tion |§| we explore different models of mass function. The 
models considered are by Press & Schechter (PS; 1974), by 
Sheth & Tormen (ST; 1999), and by Lee & Shandarin (LS; 
1999), the last two of which incorporate non-spherical col- 
lapse. In section []|, we then derive the probability distribu- 
tion function of cluster formation redshift for different mod- 
els of mass function. The formalism employed is by Lacey & 
Cole (1993, 1994), but modified to incorporate non-spherical 
collapse. In section [|, we describe how the virial mass of dif- 
ferent formation redshift is related to its temperature, and 
lay down the formalism that relates as to the cluster abun- 
dance. In section ^[ we first calculate the as normalized to 
the observed cluster abundance. We present the estimates 
of the as based on different models of mass function, and 
compare these results with other work in the literature. The 
uncertainties of our results are investigated, with an explicit 
form for the dependence on various error sources. We fur- 
ther show how to derive /3 from the information we have, 
and compare our results with recent observations. Finally, 
in section [?], we give a brief conclusion. 



2 MATTER PERTURBATIONS AND POWER 
SPECTRUM 

The standard deviation of matter density perturbations at 
a smoothing scale R is related to the mass power spectrum 
V{k) as 



a 2 (R) 



ilk 

\ w (kR)\ 2 S(k)^, 



(1) 



where w(x) is the Fourier transform of an unit top-hat spher- 
ical window: 



w{x) = 
and 
S(k) = 



3 (sin x 



2ir 2 



V(k), 



(2) 



(3) 



which is dimensionless and thus independent of the units of 
k. We shall use fe _1 Mpc and /iMpc -1 as the units of R and 
k respectively, where h is the present Hubble parameter Ho 
in units of 100 km s~ Mpc -1 . A subscript '0' will denote 
a quantity evaluated at the present epoch. Theoretically, 
the matter power spectrum Vik) in adiabatic inflationary 
models can be expressed as 



V(k) <xk n T 2 (k), 



(4) 



where n is the 'spectral index', which specifies the scale de- 
pendence of the initial perturbations, and T(k) is the 'trans- 
fer function', which transfers the initial perturbations to the 
present ep och. An accurate and a nalytically motivated form 
of T(k) is (|Eisenstein fc Hu 1998|) 



T(k) 



1 + 



[14.2 + 731/(1 + 62.5g)]g 2 
ln(2e + l.Sq) 
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r' 


X 2 /DOF (CL) 


HTP 


(o.9il l :g!) 


(o.i8l8:3) 


160+ - 085 
u.±ou_ . 051 


15.4/19 (70%) 


PD 


o.99t8:S 


r> oo-t-0.55 
u ' zo -0.16 


229+ - 042 


6.95/9 (64%) 


HTP+PD 


0.84lg;»? 


27+ - 42 


220+ 036 


24.6/30 (74%) 



Table 1. Best fits of different LSS data sets. The errors are at 95 
per cent confidence level. The bracketed numbers were obtained 
using a prior 2 > n > 0. See text for more details. 



where q = k/T (q is originally defined in units normalized 
to the inverse of the horizon size at t he epoch of radi ation- 
matter-energy-density equality), and (Sugiyama 1995) 



T = f2 m0 /iexp(-f2 B o - rT^) 

"mO 



(6) 



(•») 



Here Q m o and JIbo are the present matter and baryon energy 
densities respectively. For a fixed n, F determines the loca- 
tion of the broad peak in V(k), and thus the name 'shape 
parameter'. In the analysis of later sections, we shall use 
Obo = 0.05, and then investigate the dependence of our fi- 
nal results on this. 

From the theoretical modeling (^), we see that the 
shape of V(k) has two essential parameters, n and F, i.e., 
V(k) = V{k\n, r). Thus a comparison between this model- 
ing and observations will give us some estimates of n and T. 
The observations we considered are the linear mass power 
spectrum "P(pd)(&) by Peacock and Dodds (1994; hereafter 
PD), which was compiled from surveys of different classes 
of galaxies and galaxy clusters, and the decorrelated linear 
IRAS galaxy-galaxy power spectrum Vi(uTP)(k) by Hamil- 
ton, Tegmark, and Padmanabhan (2000; hereafter HTP), 
which was based on the IRAS Point Source Catalogue Red- 
shift Survey (PSCz; Saunders et al. 2000). As advised by 
the authors, to avoid the effect of non-linear evolution on 
smaller scales (larger k), we used only the first twelve data 
points of PD (k < 0.2/iMpc -1 ), and the first twenty two 
data points of HTP (k < 0.3/iMpc" 1 ). They are shown in 
figure |l| as squares and dots respectively (with associated er- 
ror bars). Using the maximum-likelihood method based on 
a x 2 analysis, the results of the best estimates are shown in 
table ^ and figure ^. 

In table oj, all the errors are at 95 per cent confidence 
level. The last column gives the \ 2 values of the best fits, 
the degrees of freedom (DOF) in the analysis, and their cor- 
responding confidence levels (CL). In obtaining the results 
for the HTP (the bracketed numbers), we have used a top- 
hat prior 2 > n > (see later for reasons). In figure the 
best fits of V(k) for PD and HTP are shown as the dot- 
ted and dashed lines respectively. When compared with the 
HTP (dots and dashed line), the more constraining feature 
in the PD data (squares) on larger scales (smaller k) tends 
to confine the peak of the best-estimated P(k) (dotted line) 
to smaller scales (larger k), and thus a larger estimate of F 
(see table [l]). 

We then combined both the HTP and PD into the same 
likelihood analysis (denoted as HTP+PD) with all data 
points equally weighted. In the combined analysis, we al- 
lowed a galaxy-to-mass bias in the HTP spectrum to vary: 
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Figure 1. Mass power spectra of different observations and their 
best fits (see text for details). 



"I (HTP) 



V- 



I(HTP) 



(A;) 



P, 



(HTP) 



(fc) 



1/2 



(7) 



while fixing the amplitude of V(po)(k). The best fit of 
HTP+PD showed (at 95 per cent confidence level) 



&I(HTP) — l-28jlo',18' 



(8) 



This value also implies the bias between the HTP IRAS 
galaxy power spectrum and the PD mass power spectrum. 
We note in figure jj] that for a direct comparison, we have 
rescaled down the amplitudes of the HTP data (dots with 
crosses) and their best fit (the dashed line) by a factor of 

^i(htp) = 1-28 ■ 

The estimates of Y and n from the HTP+PS are shown 
in table ^. The corresponding best fit of P{k) is plotted as 
the solid line in figure |l[ We see that it is much closer to the 
best fit of PD (the dotted line), and this is mainly due to 
the smaller error bars in the PD. From the last column of 
table [j], we also see that the best fits of each case (HTP, PS, 
and HTP+PS) are good, and agree well with each other. 
Based on the form (p|) of T, it is worth emphasizing that our 
estimates of F and n are consistent with the rece nt CMB 
constrai nts, which give 0.3 > Y > 0.1 and n ~ 1 (Balbi et 
al. 200(J; [Lange et al. 2000| ; [jaffe "et al. 2000[ ) 



One interesting observation in all the above likelihood 
analyses is that there is a strong degeneracy between the 
estimates of n and T: the larger the T, the smaller the n. 
Figure ^ shows this degeneracy in the likelihood plots of the 
three cases (PD, HTP, and HTP+PD). For any n within 
2 > n > 0, we found that the V of local maximum likelihood 
can be well fitted by 



0.657 exp(- 
0.912 exp(- 
0.868 exp(- 



1.38n 
1.38n 
1.38n 



0.02n b ) 

0.002n 6 ) 

0.002n 6 ) 



for HTP, 

for PD, 

for HTP+PD, 



(9) 



all within 3% error. This degeneracy is caused by the fact 
that V(k) is lack of feature except a broad peak, so that in 
producing the same shape of P(k) a larger F (which tends 
to shift the peak to the right) can compensate for a smaller 
n (which tends to shift the peak to the left). This is theoret- 
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Figure 2. The 68%, 95%, and 99% contours (inner out) of the 
maximum-likelihood analysis for the parameters T and n, based 
on the PD (dashed), HTP (dotted), and HTP+PD (shaded). The 
strong degeneracy between T and n is apparent. The square^ cir- 
cle, and cross label the best fit of each case shown in table 



ically expected as can be seen from equations (^) and (|J). 
Thus such a degeneracy allows us to write 



P(k;n,Y) ocP(fc;l,F'), 



(10) 



with the 'degenerated shape parameter' 

Y' = Y'(n,Y) = rexp(1.4n- 1.4) « 0.247r exp(1.4n), (11) 

which equals the shape parameter F when n = 1. We 
verified that this parameterization is accurate within 10% 
error (when V(k) is normalized at k — 0.2/iMpc _1 ) for 
1.4 > n > 0.5, 0.6 > T > 0.09, 0.28 > V > 0.15, and 
0.01 > fc/^Mpc" 1 ) > 0.3 (the range of linear scales probed 
by observations). Using this parameterization, we obtained 
the best fits of F' shown in table |l[ We note that with the 
degeneracy ( |l"l| ) the best fits of V' are consistent with those 
of the n and Y. When comparing the best fits of F' in HTP, 
PD, and HTP+PD, we also note that the larger values of Y' 
in the HTP+PD and PD analyses reflect the previously ob- 
served fact that the peak in V(k) is confined to smaller scales 
(larger k) in these two cases (see figure |l]) . In the rest of our 
analysis, we shall use the parameterization ([H^), and adopt 
the result of HTP+PD. Its errors can be well approximated 
by a log-normal distribution (at 95 per cent confidence level; 
c.f., table §) 



T' = 0.220± ;^ 



0.220 x 10 a 



(12) 



In principle, we can also use observations to constrain 
erg via equation (|l|). However, this is forbidden by the large 
uncertainty in the normalization of mass power spectrum, 
mainly resulted from the large uncertainty in the under- 
standing of the galaxy to mass bias. Hence, instead, we shall 
use the abundance of X-ray clusters to constrain erg, as is 
one of the main goals of this paper. Nevertheless, here we 
can still use Pi(htp) alone to constrain the as for the IRAS 
galaxies, namely own . A maximum-likelihood method with 
equation (113) as a prior gives an estimate 



0.78 + 0.03, 



(13) 
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Figure 3. Dependence of amplitude of mass perturbations <r on 
the smoothing scale R and the degenerated shape parameter F' . 



at 68 per cent confidence level, and crgm — 0.78 ± 0.06 at 
95 p er c ent confidence level. This result will be used in Sec- 
tion 6.3 for the discussion of redshift distortion parameter. 



We now investigate the scale dependence of the ampli- 
tude of matter perturbations. By using equation ([!]) together 
with equations (Q) and (^), we obtained a numerical fit 

g(RT) 



a(R; n = 1, r) = erg 
where as = a(R 
g(r) = r -^ T \ ijj{r) = 0.3 + 



em ' 

'<), and 



1.45 



1 + (20/r) 



Alternatively, we have 

R(<7;n=l,T)=T- 1 a— e m), 
where 

COO = Q 



t- o -0.53 

5.2s exp 



0.6 



(14) 



(15) 



(16) 



(17) 



Here g 1 means the inverse function of g, which is defined 
in equation (|l5|). These fits (eqs. [[l4| and |jii|) are accurate 
within 5% error for 1000 > RT > 0.01. We also note that 
the coupling of R and F in these fits is theoretically required 
(see eqs. jl| and |f| with the fact that q — k/T). To account 
for the dependence of cr(R) on n, we can simply replace the 
F in equation ( |l4| ) with V' (see also eq. jl^] and context), 



a(R; n, F) = a{R; n = 1, r') = a 8 



g(nr') 
e(8P) 



(18) 



The dependence of a on R and r' is shown in figure 

In the later discussion, we shall also need the red-shift 
dependence of a for a given present value in a given back- 
ground cosmology. This can be described accurately by 

<7(fim, Oa) 



<t(R,z) = a(R,z = 0)- 



(19) 



(1 + z)g(tlm0, f^Ao) ' 

where z is the red-shift, Q\ is the energy density of th e 
cosmological constant A, and (Carroll, Press & Turner 1992) 

2.5fi i i i 



g(f2 m ,OA) = 



n 



4/7 



n A + (i + fim/2)(i + n A /70) 



(20) 



which accounts for the discrepancy of the growing behav- 
ior of matter perturbations from that of a critical-density 
universe. An analytical expression for the evolution of f2 m is 



n m o(i + zf 



(21) 



(1 + z) 2 {l + Q, m0 Z - O\o) + ^A0 ' 

Throughout this paper, we shall investigate open models 
with flm < 1 and Qa = (denoted as OCDM hereafter), 
and flat models with f2 m + £7a = f2 m o + Oao = 1 (denoted 
as ACDM hereafter). 



3 THE MASS FUNCTION 

To apply constraint from the observed cluster abundance 
onto the amplitude of matter perturbations a(R), we first 
need to relate it to the number density of clusters. By defini- 
tion, the fraction of the total mass within collapsed objects 
larger than a given mass M at a red-shift z is 



T{M(R,z),z) 



fi m (> M(R,z),z) 



(22) 



ttm(z) 

Here M{R, z) is the cluster mass of a corresponding scale R 
and is defined by 

4 



M(R,z) = ^ttR 3 p(z 



(23) 



where p(z) is the average energy density of the universe 
at re d-shift z. With the numbers given in (Kolb & Turnei 



199C), we obtain 



M(R,0) = 1.162 x W 12 R 3 n rn0 h' 1 MQ, 



(24) 



where R is in /i -1 Mpc and M Q is the solar mass. Thus the 
differential number density of clusters at a mass interval dM 
about M can be derived as: 

p da 



m(M,z)dM = -~F t {p)-?-—dM, i = PS, ST, or LS,(25) 
M a aM 

where a = a(R, z), ~p = ~p~{z), /i = jtt(c), and 
cLF{M(R, z), z) 



dhia(R, z) 



(26) 



The rii in equation (^5|) is normally referred as the 'mass 
function'. In different models of mass function, the Fi(n) 
(and therefore nt) has different forms. Here we shall con- 
sider the models by Press & Schechter (1974; hereafter PS), 
by Sheth & Tormen (1999; hereafter ST), and by Lee & 
Shandarin (1999; hereafter LS): 



with [i = 

Fst(h) = 
with v = 

Fls(m) = 



/iexp i 



("*)■ 

S c /a and <5 C = 1.68647; 

0-322/| (i + riV exp (-?)' 

1.68647; and 

h) exp (^) x 

15^ 2 ^ 



(27) 



(28) 



VO.707/1, jU = Sc/a-, and 5, 
25^10 



20F 



' 5^ 
« 3 



eric(V2fi) + exp 



)+f-(-T) 
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Figure 4. The mass function rii (the differential number density 
of collapsed objects within a mass interval dM about M) of dif- 
ferent models (i =PS, ST, and LS). See equation {^) and context 
for details. 



erfc 



(49 



5v 2-7T/1 
67T 



■ exp 



(29) 



with fj, = \:j C / 'a and A3 C = 0.37. 

The n P s is based on the Press-Schechter formalism, 
which relates the mass fraction of collapsed objects whose 
mass is larger than some given threshold M, with the frac- 
tion of space in which the evolved linear density field exceeds 
some threshold 8 C . It has been extensively tested against 
N-body simulations with considerable success, but recently 
found to have steeper shape tilted towards larger mass M. 
This is mainly due to the spherical-collapse assumption 
in the PS formalism (while simulations showed that non- 
spherical collapses are important in the clustering process), 
and the fact that peaks in a linear field are poorly correlated 
with the final spatial locations of dark halos (for a review, 
see Lee & Shandarin 1999 and references therein). The ris T 
by Sheth & Tormen (1999) was first obtained as a numerical 
fit to simulations. Sheth, Mo & Tormen (1999) then showed 
that it could be associated with a model of ellipsoidal, rather 
than spherical collapse. The n L s by Lee & Shandarin (1999) 
is a full analytical alternative of n PS . It incorporates non- 
spherical dynamics and assumes that dark halos form in the 
local maxima of the smallest eigenvalue of the deformation 
tensor A3 C . It has been shown to be i n better agreement wit h 
simulations than the PS formalism (Lee & Shandarin 1999). 

Figure [| shows the mass functions m of these three 
models. It indicates that the PS mass function (dotted line) 
predicts fewer massive clusters and more light clusters than 
the ST and LS mass functions (dashed and solid lines re- 
spectively). In the literature, it has been shown that the PS 
formalism with a smaller S c can produce more massive clus- 
ters and thus resemble the simulations better for large M, 



but still su ffers from the overestim ate of light objects ( Tor- 
men 19|98|; Lee & Shandarin 1999). To show this effect, in 



figure tfwe plotted another PS mass function with S c = 1.5 
(dot-dashed line) for comparison. In addition, we should also 
note that the S c is fixed by the fit to the global mass func- 
tion, and therefore not a free parameter essentially. 



Finally, substitutingequation (^8|) together with ( |15| ) 
into (EE) (see also eq. ]is|), we obtain the present cluster 
abundance as 



rn(M,0) =Fi{n) 



3M 2 



1.45(7?F) C 



[(i?r)0-35 + 2.85] 2 J 



^1(30) 



We note that the erg dependence of m (M, 0) here has now 
been implicitly embedded into the u(<r) of Fi(fi), as well as 
the M and R through equations (hq) and (Eq). 



4 FORMATION REDSHIFT OF CLUSTERS 

While the cluster abundance in our theoretical formalism 
(eq. |3(j]) is a function of the cluster mass, the observed 
abundance of X-ray clusters is normally in the form of a 
function of the X-ray temperature. Therefore, we need to 
relate the X-ray temperature T of a cluster with its virial 
mass M v . Since clusters of the same mass may have formed 
at different red-shifts resulting in different present temper- 
atures, we first need to decompose the present number den- 
sity of clusters of given mass into contributions from dif- 
ferent formation time, and then relate the virial mass of a 
formation time with the present temperature. That is, we 
first need to find out how much of the present abundance 
rii(M, z = 0) of given mass M was formed at a given redshift 
z, and then associate this abundance with the temperature 
that corresponds to the given M and z. In this section, we 
shall discuss the first step, while leaving the second step to 
the next section. 

To achieve the first step, we need to know the probabil- 
ity that a cluster of given present mass was formed at the 
given redshift. Lacey and Cole constructed a merging history 
for dark matter halos based on the excursion set approach 
and obtained an analytical expression for the probability 
that a galaxy cluster with present virial mass M would have 
formed at the given red-shift z (Lacey & Cole 1993, 1994). 
Here we shall generalize their formalism in the following way 
to incorporate different models of mass function. At first, the 
probability that the formation time of a halo of mass Mo at 
zo was earlier than z equals the probability that it had a par- 
ent of mass M > /Mo at redshift z, where / is the fraction 
of the cluster mass assembled by redshift z. This probability 
can be easily obtained by a halo counting argument (modi- 
fied from Lacey & Cole 1993): 



P(M > fM ,z\Mo,z ) 



So 



M(S) 2(5 - So) 



dS, 



(31) 



where S = S(M) = a 2 (M,z = 0), So = 5(M ) , S> = 
5(/M ), and n = 8 c [a(M,0)/a{M,z) - 1]/V5 - 5 (the S c 
here should be replaced with A3 C for the LS formalism, and 
the same applies to the later appearance of S c in this sec- 
tion). Therefore the probability that a galaxy cluster with 
present virial mass M was formed at a given z can be ob- 
tained as: 



P*{z) 



dP(M > fM ,z\Mo,z ) 
dz 



With the change of variables 



S-Sq 
Sf - So ' 



S c [a(M,0)/a(M,z) - 1] 
\/Sf - So 



(32) 



(33) 
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equation ( |32[ ) can be rewritten as 
duo 



Pz{z) =p w (w) 



dz 



where 



p u (uj) = 



and 



dP 

duj 



Ro 
R{o) 



dFi(n) ds 



/ u= J >, <7 = «7(M,z = 0) = [a(S/ - So) + So] 



1/2 



(34) 



(35) 



(36) 



Here Ro can be obtained from equation ( J24| ) by taking 
M = Mo, and R(a) is given either by equation ( |l6| ) for 
inflationary CDM models, or by i? oc <j- 2 /(' ls +3) for power- 
law-spectrum models with P(k) oc k n " . In the latter case, we 
have [R /R{o)f = [l + fl (/-(-+3)/3 _ x j] 3/(™ s +3)^ Under 
certain conditions, equation (|35|) can be analytically inte- 
grated to yield an explicit expression of p w (w). For example, 
for power-law-spectrum models with n B — 0, we have simply 
[Ro/Ria)] 3 = 1 + s(/ _1 - 1). Thus we obtain for PS 



(37) 



p w {u) = Mf^ ~ l)erfc (75 ) - 
v/ft/" 1 - 2) exp 

and the result for ST can be also straightforwardly obtained. 

The form of equation (^H|) is general and thus suitable 
for any given model of mass function. Once given the mass 
function, which is specified by Fi(fi) (see eq. [p5[), we can 
first use equation (^) to obtain p w , and then equation ( ]34| ) 
to get p z (z). From equations (p3|), (p4|), and (|l9|), we first 
note that for fi m o = I and A = 0, w(z) is proportional to z, 
so that the conversion from p^(uj) to p z (z) requires only a 
linear rescale in both w and Even for other values of 

JlmO and A, which result in non-linear conversion from uj to z 
and thus from p^(uj) to p z (z), the relative ratios between the 
p z (z) of the PS, ST and LS formalisms will remain the same 
as in the p u (u>), i.e. the conversion from Pui(uj) to Pz(z) is 
independent of the choice of the mass function. This allows 
us to directly compare the Pu(w) of PS, ST and LS, and 
keep the conclusions still available for p z (z). Figure |^ shows 
the p w (u) in the PS (dotted lines), ST (dashed lines), and 
LS (solid lines) formalisms, for power-law-spectrum models 
with n a = —2, — 1, and 1. For a direct comparison, we have 
rescaled the p u (tii) of ST with A3 C /<5 C , and its uj with 5 C /A3 C , 
due to the fact that u> oc S c , X^c (see eq. [^]). We have also 
verified that for the inflationary CDM models, the p w (w) 
(and thus p z (z)) of a chosen mass function approximately 
interpolates between those based on the pow-law-spectrum 
models. We should also notice that the PS results here (the 
dotted lines) are identical to those presented in Lacey & 
Cole (1993). 

In comparing the results of PS, ST and LS in figure || 
we see that the p w (o>) of non-spherical-collapse models (ST 
and LS) have larger tails at high u), allowing clusters to form 
at earlier time. This is consistent with the argument that in 
non-spherical-collapse models, halos are considered formed 
as long as one of three major axes has collapsed, leading to 
the earlier formation of clusters. This earlier formation of 
clusters consumes more over-density regions at early time 
than in the PS formalism, so as to cause the delay of the 
active formation epoch, i.e. the maximum of Pui(uj) being 
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0.7 
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0.4 
0.3 
0.2 
0.1 


-0.1 
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ST 
LS 




0.5 



1.5 2 2.5 3 3.5 4 
CO 



Figure 5. The probability distribution Vi.Auj) of halo formation 
epochs u)(z) calculated from equation (p5|) of different models, 
for power-law spectra with n s = —2, —1, and 1 (from bottom 
to top at uj = 1). For a direct comparison, we have rescaled the 
Puj(uj) of ST with A3 C /<5 C , and its w with <5 C /A3 C , due to the fact 
that w on S c , \s c . Here we have used / = 0.5. We also note that 
lo(z) oc z when O m o = 1 and A = 0. 



located at a lower uj (see the solid and dashed lines in fig- 
ure ||). We also note that the dramatic increase of Pw(o;) 
towards uj — in the ST model is nearly divergent, and this 
is caused by the mathematical form of Fst- 

To implement all the analysis properly, we shall use the 
numerical result of equation ( |3q ) based on the inflationary 
models (invoking eq. flti)). The value of / we shall u se is 



(Viana & Liddle 1999; Navarro, Frenk & White 1995) (at 



95% confidence interval) 
/ = 0.75 ±0.15. 



(38) 



Thus the present number density of clusters of present mass 
M that were formed at red-shift z can now be obtained as 
rii(M, 0)pz(z)dMdz. 



5 MASS-TEMPERATURE RELATION AND 
X-RAY CLUSTER ABUNDANCE 

We then need to associate the 

abundance rii(M, 0)p z (z)dMdz with the temperature that 
corresponds the given M and z. This requires the knowl- 
edge about the relation between the virial mass M = M v of 
a formation red-shift z and its present temperature T. Here 
we shall use the result of Viana & Liddle (1996, 1999), which 
was calibrated to the hydrodynamical A^-body simulations 
of White et al. (1993b), and also shown to agree well with 
the result by Bryan & Norman (1998). Their result of the 
M-T relation is (modified from Viana & Liddle 1999): 



M v = 



(1.23 ±0.33) 



x 10 



1/2 



2(2~^)(4-^) 2 k B T 
64-56??+24i; 2 -7?7 :1 +i; 4 7.5kcV 



3/2 



f 1.67 - 
[l+zt ■ 



(39) 
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where the error is 1-sigma, z t is the turnaround red-shift, 
fi mt = n m (z t ), and 



32 



S 'A"mt 



b(Q) 



9tt 2 ^(l + zt) 3 ' 

0.76 - 0.250 (OCDM), 
0.73 - 0.23n (ACDM). 



(40) 
(41) 



For a given red-shift z = z c of cluster collapse, the 
turnaround red-shift z t is easily obtained using the fact that 
2i(zt) = t(z c ), where 



t 2 

3 



t(z) 



2(i-n m0 ) 3 / 2 
— cosh" 

■In 



-3/2 
2(1- 



for O m0 = 1, £Iao = 0, 



£ma±2)] for OCDM, 



:«2 



— 



\ 2fi m Q -ff3 m Q 

nV/ 2 +[n A o+ri m0 (i+z) 3 ] 1 / 2 
n^ 2 (i+,)3/ 2 



(42) 



for ACDM. 



We note that although other work of the M-T relation may 
suggest a different normalization from the one used in equa- 
tion ([39|), the error quoted here should have reasonably in- 
cluded the possible deviations. For example, the normaliza- 
tion used by Pierpaoli, Scott & White (2000) is about 15% 
lower than here, and this is well within our 1-sigma, which 
is about 27% of the central value. In addition, we shall also 
investigate the dependence of our final result on this, so that 
one can easily extrapolate our final result for different M-T 
normalization (see section 5.2). 



Putting all the above results together, we can now es- 
timate the present number density of galaxy clusters that 
were formed at a given red-shift z with a mean X-ray tem- 
perature k B T: 



hi(T, z)d(k B T)dz = m(M, 0)p z (z) 



3M 
2k B T 



d(k B T)dz, 



(43) 



where we have used the fact that dM/d(k B T) = 3M /2k B T. 
By integrating this over the temperature T and the forma- 
tion redshift z = z c , we obtain the theoretically predicted 
abundance of X-ray clusters observed at z b s with a temper- 
ature greater than a given threshold T t h'- 



obs 



n 1 (T,z)d(k B T)dz. 



(44) 



k B T th 



We note that the z b s on the left hand side means completely 
differently from the z on the right hand side. The former 
means the red-shift of the observed clusters, while the latter 
means the formation red-shift of those clusters relative to 
Zobs- Thus we need to integrate the above equation as if we 
were placed at z b s , i.e. we need to first shift all relevant 
physical conditions today (f2 m o, Hao, etc.) to the epoch z b s 
(fi m (z bs), ^A(zobs), etc.), and then integrate equation ( |44| ) 
as if z obs = 0. 

A comparison between the observed cluster abundance 
and the above theoretical prediction will give us an estimate 
of erg. The observed cluster abundance we sha ll use for this 
comparison is that given by Viana fc Lid dle ( |l999 ), based 
on the dataset in Henry & Arnaud (1991), and updated by 
Henry ( 2000 ). It is an abundance at z bs = 0.05 with X-ray 
temperature exceeding 6.2keV: 



AU a (> 6.2keV,0.05) 



1.53 x 10" 7±0 16 /i 3 Mpc" 3 . 



(45) 
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Figure 6. The cluster-abundance-normalized ag based on differ- 
ent mass functions. The upper and lower panels show the results 
in the OCDM and ACDM models respectively. 



The uncertainty in (^) is the 1-sigma interval, and has taken 
into account the effect of temperature measurement errors. 
The reasons for concentrating on galaxy clusters with tem- 
perature larger tha n 6.2k eV has been extensively discussed 
by Viana & Liddle $99$ ). 



6 RESULTS AND DISCUSSION 

6.1 The erg 

By comparing the observed cluster abundance ( plB] ) with the 
theoretical prediction (^), we obtained the values of erg 
based on different models of mass function. Figure ^| shows 
the results as functions of f2 m o and SIao- These results can 
be fitted by 

a s{l) {n m0 ,tt A0 ) = cifCo, i =PS, ST, LS, and ST+LS, (46) 
where 



a = a(f2 m0 , Oao) = —0.3 



0.17ns, 



o.i3n AO , 



(47) 



and the values of ci and C2 are given in table ^. These fits are 
accurate within 2% error for 1 > n m o > 0.1. We notice that 
these results of erg have only simple dependence on Qa (the 
last term of eq. ^7|), and this dependence is independent 
on the choice of mass function. It is clear that for a given 
n m o, the cosmological constant tends to increase the nor- 
malization of erg . Also given in equation (^) and table |^ is 
the mean of the ST and LS results, labeled as 'ST+LS' (see 
next section for more details). It is plotted as the thick solid 
line in figure 0, and we shall refer to this as the combined 
result of non-spherical-collapse models. 
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i 


PS 


ST 


LS 


ST+LS 


Cl 


0.54 


0.50 


0.455 


0.477 


C2 


0.45 


0.37 


0.31 


0.34 



Table 2. Values in the fit ( |46| ) of <jg, based on different models 
of mass function. 



Figure |fj indicates that the cr g (ST) and (Tg(LS) are system- 
atically lower than <r g ( P s). On average within the ranges of 
QmO and SIao probed here, we found that (Ts(rt) is 7% and 
°"8(ls) is 15% smaller than <r g (ps). In figure M, we also see 
that the combined non-spherical-collapse results o" g (ST+LS) 
(thick solid lines) are on average 11% lower than cr g (pg) 
(thick dashed lines). This implies that the inclusion of non- 
spherical collapse considerably reduces the amplitude of 
mass perturbations required for resembling the observed 
cluster abundance. This can be understood as the conse- 
quence of the higher abundance of massive clusters in the 
non-spherical-collapse models. In figure ^j, we saw that the 
non-spherical-collapse models (ST and ST) predict more 
massive clusters. Adding the fact that clusters are relatively 
rare objects that exist only with high mass, we know that the 
amplitude of matter perturbations required to resemble the 
observed cluster abundance will be lower in non-spherical- 
collapse models, resulting in the lower <r g as seen in figures ^ 
and ^ 

In figure [?], we have also plotted other results in the lit- 
erature for comparison. They are results by Viana & Liddle 
(1999; VL, thin solid lines), by Eke, Cole & Frenk (1996; 
ECF, thin dashed lines), by Wang & Steinhardt (1998; WS, 
dot-dashed lines), by Borgani et al. (1999; B, dotted lines), 
and by Pierpaoli, Scott & White (2000; PSW, circles in the 
lower panel). All these results were based on the PS formal- 
ism except for the last one (PSW), which was based on the 
ST mass function. We see that for fi m o > 0.2, our cr g (ps) 
is consistent well within 10% error with all the PS-based 
results except for the Borgani et al. (1999) in open mod- 
els. In particular, our cr g (ps) is consistent with that of VL 
within 5% error. Since our PS framework is parallel to that 
of VL, we verified that this small deviation is partly caused 
by the fact that we used the inflationary CDM model to 
calculate p z (z) (the probability that a cluster of given mass 
was formed at z), while VL used thep z (z) based on a power- 
law-spectrum model with n s = (i.e., our eq. J§). On the 
other hand, for non-spherical-collapse models, our er g (ST+LS) 
is lower than that of PSW. This discrepancy is mainly caused 
by the fact that we integrated the abundance over the for- 
mation shift while the PSW did not. We verified that the 
ignorance of the integration over the formation shift will re- 
sult in a significant overestimate of <7 8 , especially for low 
fi m o cosmologies. This is because clusters of the same mass 
but formed in the past will have higher temperature and 
thus more observable than those formed today, so that the 
required <j g to resemble the observation in models allowing 
clusters to form in the past will be lower. We also note that 
our normalization for the M-T relation is higher than that 
of the PSW. If we adopt their normalization, our o- 8 (st+ls) 
will be even lowered by about 5%. Overall, it is clear that 
our ct 8 (st+ls) is lower than any previous results, and this is 




° i> 6 

m0 



Figure 7. The cluster-abundance-normalized ag based on 
ST+LS (thick solid lines) and PS (thick dashed lines), compared 
with results of Viana k. Liddle (1999; VL, thin solid), Eke, Cole & 
Frenk (1996; ECF, thin dashed), Wang k, Steinhardt (1998; WS, 
dot-dashed), Borgani et al. (1999; B, dotted), and Pierpaoli, Scott 
& White (2000; PSW, circles in the lower panel), in the OCDM 
(upper panel) and ACDM (lower panel) models. The shaded ar- 
eas are the 1-sigma regions of the ST+LS result. The crosses are 
the observational constraints derived by Henry (2000). 



mainly attributed to the inclusion of non-spherical collapse 
in our formalism. 

It should be also noticed that in Henry (2000), the max- 
imum likelihood fits of cluster evolution models to the ob- 
served cluster X-ray temperatures gave f2 m o = 0.49 ± 0.12 
and cr g = 0.72 ± 0.1 for OCDM, and « m0 = 0.44 ± 0.12 and 
<r g = 0.77 ± 0.15 for ACDM, all at 68 per cent confidence 
level. These results are plotted as the crosses in figure ^. 
They are consistent with all the results presented here. 



6.2 Uncertainties 

To know how stringent our results of a$ are, we investigate 
the dependence of <r g on M v , f,N= N ohs (> 6.2keV,0.05) 
and F', because they carry uncertainties in our analysis 
pipeline. Empirically we found 



<7 8 (r',/,M v ,iV) 

°"8 (r' ( 0) , /(o) , M v(0 ) , -/V(o) ) 



v(0) 



P_ 

r' 

(0) 
PM 



7(0) / 

N + N 
N(p) + N 



pi- 



rn 



where the subscript (0) denotes the central values (see 
eqs. Q, 01, H, and 0) that were used to obtain the 



main result 
p r , = 0.14 



and 



■ 0.13O23 +0.1J&), 



(49) 
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pf 

Pm 



p N = 0.12n n 



1.95-1.9fCo -0.25«ao, 
0.41 - 0.06^0+0.08^0, 

+ 0.05Q 2 AO - 5p S 0.02in: 



N = 5 x 10~ 8 /i 3 Mpc~ 3 , 



(50) 
(51) 
(52) 
(53) 



with Sps equal unity for the PS, and zero otherwise. Equa- 
tion is accurate within 3%, 3%, 2%, and 0.5% errors 
for 



(54) 

(55) 
(56) 

(57) 



1.8 > 


M v 


> 0.2, 


Afv(0) 


0.9 > 


/ 


> 0.53, 


3 > 


N 


> 0.3, 


0.28 > 


r 


> 0.17. 



From equation (|48|), we first note that within the parameter 
range probed here (1 > Qmo > 0.1 and Qao = 1 — fimo or 
0), the as is always a monotonically increasing function of 
all A/ v , N, f and F'. Second, it can be also inferred that 
within the 68% confidence level of M v , N, f and F' (see 
eqs. Q, Q, Q, and @), the resulting as within the 
parameter range probed here can vary up to 14%, 6%, 5%, 
and 2% respectively, depending on the background cosmol- 
ogy. Fhus we see that an improvement in reducing the un- 
certainty in M v , the normalization of the M-T relation, can 
most efficiently reduce the uncertainty in the resulting as- To 
the contrary, a more accurate estimation of the degenerated 
shape parameter F' will not change the resulting erg much. 
This observation helps relax the worry about the accuracy 
in using the parameterization (^), and the approximation 
([]), which we have used to model baryonic effects. 

To investigate the overall uncertainties in as (that is, 
to include the uncertainties from all M v , f, N and T' si- 
multaneously), we implemented Monte Carlo simulations, 
with the distributions of the uncertainties in V' and A r Q b s (> 
6.2keV, 0.05) as log- normal, and those in M v and / as Gaus- 
sian (see eqs. (T§, @, |§, and ||). We found that the 
distribution of the resulting as is very much Gaussian, with 
a mean given by equation (^) and table [| and a standard 
deviation (68% confidence level) 



6i(f2m0, f^Ao) = 0.1f2 m0 ' a S (i) (fimO, f^Ao), 



i = PS, ST, LS.(58) 



We have verified that this result is very weekly dependent 
on Qao and the choice of the mass function (PS, ST, or LS), 
and these two aspects together contribute only a deviation 
of less than 3% from the above result. 

To obtain a single result for the non-spherical-collapse 
models (ST and FS), we linearly combine the probability 
distribution functions of ag(ST) and cr g (xg) obtained from the 
above Monte Carlo simulations. The resulting distribution 
of as is very Gaussian. It has a mean given by equation ( |46| ) 
and table H (the ST+FS result), and a standard deviation 

3 ' 12 °"8(ST+LS)(^m0, ^Ao) 



£ST+Ls(fimO, ^Ao) = 0.11£l m0 ' O"g(ST+LS)(^m0, ^Ao)- (59) 

Figure ^ shows this ST+FS result as the thick solid lines 
(the mean; eq. ]46| and tab. ^ with shaded areas (the 1- 
sigma regions; eq. |59|). Also plotted for comparison is the 
PS result (the thick dashed lines). It is clear that the result 
based on non-spherical-collapse models is about one sigma 
lower than the PS result, due to the origin already argued 



in section 3.1 



verified that for 0.05 > SIbo > 0, the normalization of as 
was affected by less than 1%. We have also tested that the 
difference between the as based on the p z (z) of the stan- 
dard CDM model and of the power-law-spectrum model 
with n B = is less than 3% in all cases. In addition, we 
explore the effect of the integration over the formation red- 
shift p z (z) based on different models of mass function, i.e. 
while keeping the mass function the same (PS, ST or FS), 
we use different p z {i)(z) (i =PS, ST or FS) to calculate 
as- It is found that for a given model of mass function, 
the erg is higher when using p z (LS){ z ) than using p z (ps){z), 
and is even higher when using p z (ST){ z )- Adding the fact 
that p z (ST){ z ) and P z (ls)(z) have maxima located at lower 
z than p z (ps)(z) (see figure ||), the above result is consis- 
tent with the previously argued and verified fact that the 
required <rg to reproduce the observed cluster abundance 
will be lower if clusters are actively formed earlier (and thus 
have higher temperature and are more observable). Never- 
theless, the difference in the erg using different p z (z) but the 
same mass function is always less than 3%. Therefore, we 
conclude that the previously observed lower normalization 
of as in the non-spherical-collapse models is indeed mainly 
caused by the form of the mass function itself, rather than 
by its form of the formation-redshift probability, which in 
fact has opposite effect. 

6.3 Redshift distortion parameter 

Another important aspect in the study of large-scale struc- 
ture is the so-called 'redshift distortion parameter'. It quan- 
tifies the confusion between the Hubble fl ow and the pe culiar 
velocities, and is analytically defined as (Kaiser 1987) 



Pj (flmO, ^Ao) = 



/(f2 m 0, Oao) 



= /(HmO, f^Ao) 



0~8 



(60) 



where bj = as(j)/as is the 7-type galaxy to mass bias 



(j =IRAS, optical, etc.), and ( Fahav et al. 1991 ) 
/(fi m o, ^ao) 



n 0.6 . ^A0 m fi m o. 

" m0 + — (1 + — ) 



(61) 



which is the rate of growth of matter perturbations at the 
present epoch. Here we shall estimate (3i, with the subscript 
T indicating the IRAS galaxies. To this end, we can substi- 
tute the result (f|]) (which is calibrated from cluster abun- 
dance) and the result (|l^) (which is obtained from the IRAS 
PSCz survey) into equation (|60|). By combining their likeli- 
hoods together, we obtained the maximum-likelihood result, 
which can be fitted as 



d 2 -0.16(fi m0 +fi A0 ) 



=PS, ST+FS, (62) 



where the values of d\ and d,2 are given in table |3|. This 
fit is accurate within 2% error. The errors at 68 per cent 
confidence level was also found as 



£/3(i)(^m0, S7ao) = 0.176f2 mQ ft(i)(ftmO, f^Ao). 



(63) 



In our entire analysis, we have used f2p,o = 0.05, but 



This fit is accurate within 1% error. These results are plotted 
in figure ||. 

From figure^], we first note that the /3i(st+ls) is system- 
atically lower than /?i(ps). This is a direct consequence of the 
fact that /3i oc as, which is lower in non-spherical-collapse 
models as previously observed. We also see that the inclu- 
sion of the cosmological constant tends to increase f3i. In the 
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i 


PS 


ST+LS 


di 


0.693 


0.613 


d 2 


0.26 


0.24 



Table 3. Values in the fit (g2j) of Pih\, based on different models 
of mass function. 
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Figure 8. Linear redshift distortion parameter /3i, estimated from 
the combined maximum-likelihood analysis of &s(i) and erg. The 
solid lines are based on non-spherical-collapse models (ST+LS), 
with the shaded areas representing the regions within 68 per cent 
confidence level. The dotted lines are based on the PS formalism. 



different models of mass function. We found that clusters of 
the same mass can form earlier but have a later epoch of 
active formation in the non-spherical-collapse models (ST 
and LS) than in the PS formalism. This is consistent with 
the observation from numerical simulations. Based on dif- 
ferent models of mass function and their associated proba- 
bility distributions of formation redshift, we then used the 
observed cluster abundance to estimate the amplitude of 
matter density perturbations on the scale of 8/i _1 Mpc, as, 
and the redshift distortion parameter for IRAS galaxies, Pi, 
in both OCDM and ACDM cosmologies. The <rg and Pi re- 
sulted from non-spherical-collapse models are systematically 
lower than those based on the PS formalism. We showed that 
this is mainly owing to the larger mass function at the high 
mass end in the non-spherical-collapse models. The origins 
of the uncertainties in our final results were also investi- 
gated separately, and we found that the main contribution 
is from the uncertainty in the normalization of the virial 
mass-temperature relation. Therefore we expect that fur- 
ther improvement in the study of this normalization will 
provide us with more stringent constraint on both erg and 
Pi. In addition, since we saw significant corrections in the 
resulting erg and Pi when switching from the conventional 
Press-Schechter formalism to the more recent non-spherical- 
collapse models, we urge the use of these models in all rele- 
vant studies, especially when we are entering the regime of 
precision cosmology. 
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past, pi has been measured by many authors using different 
observations (for a review, see Strauss & Willick 1995) . The 
combined result of these observations has a central value of 
about Pi ss 0.7, but with a large 68 per cent confidence re- 
gion of about ±0.4. In recent years, observations seem to 



favor a lower value. For example, Pi = 0.41 



+0.13 



by Hamil- 



ton, Tegmark & Padmanabhan (2000), and fa = 0.39T0.12 
by Taylor et al. (2000). If these measures are correct, then 
the PS-based /3r(ps) (dotted lines in fig. ^) will be in trouble 
except in the very low f2 m o regime, and the non-spherical- 
collaps Q models will serve to relax this situation (solid linos 
in fig. 



7 CONCLUSION 

In this paper, we first used the observed linear mass power 
spectra to estimate the spectral index n, the shape param- 
eter r, the degenerated shape parameter F', and the IRAS 
<t 8 (i) (see table [l], eqs. jy], |13)). We then derived the prob- 
ability distribution function of cluster formation redshift for 
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